Decoherence from a spin-chain with three-site interaction 
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We investigate the time evolution of quantum discord and entanglement for two-qubit coupled to a spin 
chain with three-site interaction in the weak-coupling region. If the quantum system evolves from a Bell state, 
quantum correlations decay to zero in a very short time at the critical point of the environment. We found 
there exist some special interval of the three-site coupling strength in which the decay of quantum discord and 
entanglement can be delayed. When the qubits are initially prepared in a Bell diagonal state, the decay of 
entanglement is also delayed in the special interval, but the decay of quantum discord is enhanced. Besides, the 
sudden transition between classical and quantum decoherence is observed, we found the transition time can be 
lengthened at the special range of three-site interaction and shorten by degree of anisotropy. 

PACS numbers: 



I. INTRODUCTION 

Entanglement, which is an important characteristic trait of 
quantum mechanics, is playing a key role in quantum informa- 
tion processing [1]. It is a type of quantum correlation arising 
from the superposition principle of quantum mechanics and 
has a wide range of applications on quantum teleportation 
quantum cryptography [3], and universal quantum computing 
H. In fact, it is only one aspect of quantum correlation and 
there exist other nonclassical correlations |5-9]. Recent stud- 
ies have indicated non-entanglement quantum correlation can 
be responsible for the quantum computational speed of deter- 
ministic quantum computation with one qubit |10]. Quantum 
discord was introduced by Olliver and Zurek to quantify the 
nonclassical correlation of bipartite system and it is defined 
as the difference between two quantum versions of classical 
mutual information |5|] . 

Quantum discord has been studied widely, such as its dy- 
namical property and indicating quantum phase transition 
lfTol - [T6ll . It has been shown thermal quantum discord in 
Heisenberg models is different from thermal entanglement in 
many unexpected ways, and the results point out that it may 
be useful in the experimental detection of critical points for 
quantum phase transition 1 1 1 ] . Quantum discord could be 
considered as a figure of merit for the quantum advantage in 
some computational models with little or zero entanglement 
lfTol[T2ll . T. Werlang has pointed out that quantum discord is 
more robust than entanglement against decoherence and con- 
cluded quantum algorithms based on quantum discord can be 
more robust than those based on entanglement [lit]. Ref. jji}] 
show that quantum discord may reveal more general informa- 
tion than quantum entanglement for characterizing the envi- 
ronment's quantum phase transition. 

For a realiatic quantum syatem, it is inevitable to interact 
with the environment and the interactions between the quan- 
tum system and its environment leads to decoherence which 
leads to the destruction of quantum properties. The deco- 
herence behaviors of correlation have been investigated in 
Markovian and non-Markovian environment and the results 
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show that quantum discord vanishes only in the asymptotic 
limit in all cases under the Markovian environments |13] and 
instantaneously disappear at some time points under a non- 
Markovian environment (vf\ . Moreover, it is found either the 
classical or the quantum part of correlation may be unaffected 
by decoherence fl8Lll9ll . In addition, the dynamic behavior of 
quantum discord under spin environment has been discussed 
ll20[|2Tll . It is found that the quantum discord for the two cen- 
tral qubits can become minimized rapidly close to the critical 
point of a quantum phase transition |,2ftl . Besides, the DM in- 
teraction is considered in spin environment. The results imply 
the DM interaction can enhance the decay of quantum corre- 
lations and this effect can be strengthened by anisotropy pa- 
rameter Recently, studies of decoherence induced by 
spi n environment have received a great deal of attention |22- 
12411 and the results revealed that quantum coherence is dra- 
matically destroyed at the critical point of the quantum phase 
transtion of the external environment ll23ll . 

Previous studies have been shown that the three-site inter- 
action may play an important role in quantum entanglement 
and discord of the bipartite system ll25[ - l27ll . Thus, it is nec- 
essary to consider when environment presents three-site inter- 
action, how it affects the evolution of quantum discord of the 
central two-qubit system. In this paper, we offer a compara- 
tive study of the dynamic behaviors of quantum discord and 
entanglement for two-qubit prepared in a class of X-structure 
state under spin chain environment. We will show how the 
three-site interaction affects those process of the nonclassical 
correlation. The paper is organized as follows. We start in sec- 
tion II by introducing and diagonalizing the model, and obtain 
the reduced density matrix of the central two-qubit system. In 
section III, we analytically and numerically evaluate quantum 
discord and entanglement, then present the main results. The 
last section is devoted to the conclusions. 



II. MODEL AND ITS SOLUTION 

The total Hamiltonian for two central qubits coupled to an 
XY spin chain with three-site interaction is composed of two 
parts: the interaction Hamiltonian H; of the central two-qubit 
with the chain and the self-Hamiltonian He of the spin-chain 
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environment. They can be written as 



agonalized Hamiltonian as 
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where denotes coupHng strength, cr^^^j and cr^(a - 

X, y, z) are PauH matrices of the two central qubits and the sur- 
rounding chain, respectively. is the total number of the sites 
of the chain. The parameter A characterizes the strength of 
the transverse field, y defines the degree of anisotropy of the 
interactions in the x - y plane, and a denotes the strength of 
three-site interaction. 



The eigenstates of the operator (^cr^ 



given by 



- ^cr^) are simply 



I0l) = |00),|<^2) = 101), 103) = 110), 104) = 111), (2) 



and the Eq. (1) can be written as 



H 
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where the spectrum reads 



A7 = 2 
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(5) 



in which e^" = {[cos (M) _ ^^j' + ^2 ,i„2 



With the calculation above we can obtain the time evolution 
operator for the Hamiltonian: 
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t/(f) = j]|0;.)<0;.|®f/£''(O, 



(6) 



^=1 



where U ^ (f) = exp [-iH^'i^ denotes the projected time evo- 
lution operator. Then the reduced density matrix of two cen- 
tral qubits is obtained by tracing out the environment, 

PABit) = Tr£[p'°'(f)] 

4 

= YjP^v{t){cl)^\pABm<Pv)\(l>f.){(l>v\, (7) 

where F^y (t) = {(fslU^/' (t) t/^" (t) \ipe)- Now we assume that 
the two central qubits are initially prepared in an X-structure 
state 



where the parameter /l^ is /li(4) - A + g, /l2(3) = A + g6, and 



is obtained from by replacing A with /l^. The cen- 



tral two-qubit system and the chain environment are supposed 
to be initially uncorrected with the density operator p'°'(0) - 
Pab(O) ® \<f)EEW\^ where Pab(O) and |(/?)£ are the initial states 
of the two central qubits and the spin chain environment, re- 
spectively. The time evolution of the total system is governed 
by p'°'(f) = U(t)p'°\Q)U\t) with U{t) = sxpi-iHt). To ob- 
tain analytical expression of the time evolution operator, we 
need to diagonalize the projected Hamiltonian H'^' . With the 
help of Jordan- Wigner transformation given by 



Y](l-2clc,)(c, + cl), 
-iY](l-2clc,)(ci-c]), 



s<l 

1 - 2cjci, 



where c; and c| are spinless fermion annihilation and cre- 
ation operator respectively and employing Fourier transforms 
of the fermionic operators described by dk - -^l.icie^'^""'^^ 
with k = -M - ■ - M and M = (A^ - l)/2, following by Bo- 



goliubov transformation yiij^^ = cos -|-iiA- - ' sin -Yd_^. with 
0^" = arctan [y sin (^) / \a^ - cos (2^)]}, we obtain the di- 



Pab (0) = 



1 



(8) 



where m - 1,2,3 and / is the identity operator In order to 
expediently obtain an analytical result, we set the parameters 
c,„ to be real. According to Eqs. (7) and (8), in the basis 
{|00), |01), |10), 111)), the reduced density matrix of the two 
central qubits is given by 



n H- C3 



Pab(0 






G* 



G 
1 - C3 W 
1 -C3 

1+ C3,j 



(9) 



where G = (ci - C2)Fu(t), W ^ {c, + C2)F23(t), F^y(t) 
is the decoherence factor, and * denotes complex conju- 
gation. It is necessary to obtain the exact expression for 
the decoherence factor. Then the groundstates of pure 
spin-chain Hamiltonian H'^ and the qubit-dressed Hamilto- 
nian H'^^ can be denoted with |G)^ and |G)i^. We as- 
sume that the initial state of the environment |0)£ is the 
ground state |G)^ of the pure spin-chain Hamiltonian, |G),i = 

®jJlj(cos y)|0)i|0)_i. + i sin Y)\^)k\'^)-k- By using the transfor- 
mation \Gh = nf=i(cos0^'' + isin®lyl^y_^ -^)\Gh^ with 
0^^" = iol" -61)12, We can obtain the decoherence factor IH, 
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FIG. 1: (a) Entanglement and (b) quantum discord as a function of 
time t and three-site interaction a. Other parameters are set as 7= 1 , 
6 = Q,g= 0.05, i= 1, and W = 400. 



FIG. 2: Quantum discord (solid line) and entanglement(dashed- 
dotted line) as a function of time t for different a. Other parameters 
are set as 7= 1, (5 = 0, ^= 0.05, a= 0, and = 400. 



M 

\F^At)\ = - sin2(20^'') sin2(A,'''f) - sin2(20f ) 

X s,m\A.ff) + 2 8111(20^") sin(20f ) 
X sin(A^''f) sin(A^' f) cos(A^''f - Af f) - 4 
X sin(20^'") sin(20f ) sin2(0^'' - 0^') 
xsin2(AS)sin2(A;'''0]'/l 



(10) 



When F\ii(f) — > 0, it denotes that the quantum coherence 
between the two central qubits is strongly destroyed and un- 
dergoes strong decoherence due to the interaction with the 
chain environment. When Fi4(f) — > 1, the two-qubit system 
is slightly affected by the environment. 



III. CORRELATIONS FOR THE TWO-QUBIT 

Quantum discord is defined as the difference between the 
quantum mutual information and the classical correlation 
shared by two subsystems and expressed as: 



QipAB) = I{Pab) - C(pab) 



(11) 



where I{pab) is the total correlation and C (pab) represents 
the classical correlation of state pab- These two correlation 
can be written as: 



I(Pab) = S{pa) + S(pb) - S(pab), 



(12) 



S(Pab\\^^]) = ^jPjSiPAin") is a quantum extension of the 

classical conditional entropy, and pAin'^ is the reduced state of 

j 

partition A after the measurement 11^ performed on partition 
B, with outcome j. 

To obtain the quantum discord of pABit), we need to calcu- 
late the quantum mutual information and classical correlation. 
It is not difficult to calculate the quantum mutual information, 



IlPABiO] = 2 + ^ w„ log2 a>„. 



(14) 



where wi 



1, 



n=l 



|W|),W3 = i(l+C3 + 



H = iil-C, + \W\),C^2 = ^(1-C3 

|G|), and 0)4 - ^(1 + C3 - |G|). To calculate the classical corre- 
lation C{pab), we take the complete set of von Neumann mea- 
surement n^' = |<^*^')(0'-'*| as a local measurement performed 
on the subsystem B, where = cos0|O) + e"^ sin 011), 

= e""^sin0|O) - cos 011) are the two projectors with 
6,(p e [0, 27r]. Thus we obtain the reduced density matrices 
of subsystem A after measurement. 



^[1 



• C3 cos(26')] 



i[l -C3 cos(26i)] 



(15) 



(2) /i[l-C3Cos(20)] -r \ 

\ -r i[l +C3COS(20)] j' ^""^ 

and the probability pi - p2 - where F - j{e^"^W + 
e'^G) sin(20). Subsequently the classical correlation of Eq. 
(9) can be calculated as 



C(pab) = max{5(pA) - S {pABm^j))), 



(13) 



where Sip) - Trplogjp is the von Neumann entropy, Pa{B) - 
TrB(A)(pAB) is the reduced-density operator of partition A{B), 



C\pAB{t)] = ^ log2(l +x) + ^ logzd -X\ (17) 



where - max||c3|. 



W+IGll 



Finally, substituting Eqs. (14) 
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FIG. 3: (a) Entanglement and (b) the quantum discord as a function 
of time t for different a. Other parameters are set as y= \, S = Q, g= 
0.05, A= 0, and A' = 400. Here a = Q (solid line), a = 0.1 (dashed 
line),^ = 0.5 (dotted line) 



and (17) into Eq. (11), we can obtain the expression of quan- 
tum discord immediately. 

In order to compare the dynamical behaviors between quan- 
tum discord and entanglement, we need to calculate entangle- 
ment. Now we choose entanglement of formation to quantify 
the amount of entanglement for the two-qubit. 



1 + 



(18) 



where /(jc) - -.xlogj .x-(l -.x) log2(l -x) is the binary Shan- 
non entropy, and C - maxjO, 2A,„ax - Tr{ -\/Pab£?J) is the time- 
dependent concurrence [28^ with g - <^\®crBPAB'^\®'^B' 
^max is the the maximum eigenvalue of yfpABQ- We can 
obtain concurrence for the state given by Eq. (9), C - 
inax{0,^^,ffl^}. 



A. Evolution from Pure State 

The dynamical behavior of quantum discord and entangle- 
ment for two qubit in XY spin environment have been dis- 
cussed in 1 20]. Here we focus on the dynamic phenomena of 
the two-qubit due to the three-site interaction. Let us consider 
the two-qubit with pure initial state. The parameters are set 
as c\ - -C2 = C3 - 1, then the initial state becomes Bell 
state -;^(|00) + |11)). Substituting these values into Eqs. (11) 
and (18), we obtain the quantum discord and entanglement 
of this case, Q = \og^( l - IF^I) + ^ log^d + |f wl). 



1- ^/hl 



■l0g2(l- 



^/T^-\Fl4\^). Obviously, both of them are monotonically in- 
creasing function of variable |Fi4(f)l in the interval [0,1). 
Thus, they behave in a similar way. 

Fig. 1 shows the effects of the three-site interaction, which 
is consistent with the result in f2^, where the concuiTence is 
discussed in the system. The dynamical behaviors of quan- 
tum discord and entanglement are similar as expected and 
they both decrease monotonously against time. When t is 



FIG. 4: (a) Entanglement and (b) quantum discord as a function of 
time / and the three-site interaction a. Other parameters are set as 7= 
1, ^=0, g= 0.05, A= 1, and A' = 400. 



fixed, quantum discord and entanglement behave as parabo- 
las and have maximum value at the same symmetry axis 
a - a'. When a is fixed, quantum discord and entangle- 
ment decay to zero in a short time. Furthermore, one can 
observe that the decay of quantum discord and entanglement 
become more sharply with increasing intensity of three-site 
interaction a in the region a > a' and with decreasing a in 
a < a'. At the point a = a', the decay of quantum dis- 
cord and entanglement be delayed notably. To understand 
this effect. One may turn to the approximation of \Fi4 (t) \ 
given in Ref. l26ll . For a weak-coupling parameter g and 
a large A^, |Fi4(f)l ~ g-f^'+^^+^^V- with n = ^j^-Lij-f, 
X, ^ 25fa-VV ^(|)3^ ^ lii£WE(*)4. When we fix f, 

(/i— 1) ^ (/I— 1) i-' 



the exponent (ti +T2 + T3) is a quadratic function about a, so 
there should be an extremum point where \Fi4 (f) | reach the 
maximum value. Since quantum discord and entanglement 
are monotonically increasing function of variable \Fi4 (f) | in 
the interval [0, 1), they must reach maximum values at the ex- 
tremum point. Consequently, the decay of quantum discord 
and entanglement be delayed at a'. 

Moreover, from the expression of \Fi4 (t) |, we can see when 
/i — > 1, 1^14 (f) I will gradually decay against time. Fig. 2 
shows the influence of transverse field strength A on the de- 
cay behavior of quantum correlations. One can observe that 
quantum discord and entanglement decrease against time as 
expected before they tend to stable values, when A = 1 quan- 
tum discord and entanglement suffer sudden death at the same 
time. 

The situation of A - I have been considered in Fig. 1 and 
we found there exist special interval where the decay of the 
quantum discord and entanglement can be delayed remark- 
ably. Here we consider the case of A = 0. We plot quantum 
discord and entanglement as a function of time t with differ- 
ent three-site interactions at /I = in Fig. 3. It is found that 
the dynamic behaviors of quantum discord and entanglement 
are obviously distinct from the case of A - I. The behaviors 
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FIG. 5: Quantum discord and entanglement as a function of time t 
witfi different a, Otiier parameters are set as A= 1, (5=0.5, g= 0.05, 
A= 1, and A' = 400. Here, a = -0.8 (solid line), a = -0.5 (dotted 
line), a = (dashed line), a = 0.5 (dashed-dotted line). 

show oscillatory decay against time. From Fig. 3, one can 
find the absolute value of amplitude is very sensitive to the 
three-site interaction. In other words, the absolute value of 
amplitude is decreasing with increasing three-site interaction 
a, and its decay can be enhanced by a. However, the aver- 
age value of quantum discord and entanglement is approxi- 
mately remain unchanged. What's more, the time between 
the adjacent peak is shorten with increasing the value of a. To 
explain these features, we take a similar tactic employed in 
Fig. 1. |Fi4(f)l = where d = -^^,^^^^1^ ^l' k^, 

w = 1 -H Obviously, the exponent of |Fi4(f)l is quadratic 
sine function about a, and this finding indicate attenuation vi- 
bration of the nonclassical correlations, which is consistent 
with the numerical results shown in Fig. 3. 

B. Evolution from Mixed State 

In this subsection, we consider the case of two-qubit ini- 
tially prepared in the mixed state. We choose the parameters 
ci - I, -C2 = C3 = 0.2, so the initial state becomes pab (0) - 
O.6|<^+)(0+| + 0.4|^+)<(^+|, where = ^ (|00> -i- 111)) and 

|^-^> = ^(|01) + |10)). 

In Fig. 4, we show quantum discord and entanglement as 
a function of a and f. One can see that in the special interval 
nearby a - -0.5, the decay of entanglement can be delayed 
remarkably, while the decay of quantum discord is enhanced. 
To understand these features, we take a similar tactic em- 
ployed in Fig. 1. Q = 0.41og2(1.6)-(-0.2-i-A:i log2 A'i-A-2 log2 A'2, 
E - -/c3log2/f3-(l-A:3)log2(l-/f3), whercA-i = 0.3-1-0. 3|Fi4|, 

m , mir i 1+ Vl-[0.6|Fi4(f)|-0.4p „, . , 

K2 - 0.7h-0.3|Fi4|, /C3 - — 5 , Obviously, quan- 
tum discord is a monotonically decreasing function of variable 
\Fi4 (t) I, while entanglement is monotonically increasing at 
interval \Fu (t) | e [0.667, 1] and zero at |Fi4 (t) | G [0, 0.667). 
Since there be an extremum point where \Fi4 (f) | reach the 
maximum value. So at the extremum point, the decay of quan- 
tum discord is enhanced while the entanglement can be de- 
layed. Moreover, from the plot we find that while the entan- 
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FIG. 6: Quantum discord (solid line), classical correlation (dotted 
line), and total correlation (dashed line) as a function of time t. Other 
parameters are set as o- = 0, y= 1, 6=1, g= 0.05, A= 1, and N = 400. 
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FIG. 7: The transition time as a function of three-site interaction a. 
Other parameters are set as 7= 1,6=1, g= 0.05, A= 1. 



glement suffers sudden death after a finite time, the quantum 
discord is maintained a stable value during the whole time 
period. This imply that the nonclassical correlations may be 
more resistant to external perturbations than that of entangle- 
ment. 

The obvious difference between the behaviors of the quan- 
tum discord and entanglement can be found in Fig. 5, where 
we consider the influence of the anisotropy parameter One 
can observe that quantum discord and entanglement behaves 
in a different way. Quantum discord increase at the first sev- 
eral seconds and then decay with time, while EOF only decay 
monotonically. At the point of a = 0, the decay of entangle- 
ment can be enhanced by increasing parameter 7, while the 
increase of quantum discord is strengthened, and the decay 
is weaken. For the case of a 5^ 0, we observe that the de- 
cay of entanglement can be delayed by increasing the value 
of y. However, under the same situation, both the decay and 
increase of quantum discord are notably enhanced. What's 
more, entanglement quickly vanishes in a much shorter time 
than quantum discord does. 

We now turn to investigate the transition behavior of quan- 
tum discord which has been discussed in Ref. ll29ll . In Fig. 
6, we plot the dynamic behaviors of the quantum discord, the 
classical correlations and the mutual information against time. 
We set the parameter 6-1, which means one qubit interacts 
with the surrounding spin chain and the other is free of the 
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FIG. 8: The transition time as a function of the degree of anisotropy 
y. Other parameters are set as a= 0, (5=1, g= 0.05, X= 1, and A' = 
400. 

environment. The plot clearly shows the sharp transition from 
the classical to quantum decoherence and the transition occurs 
at f = f'. In the first f' seconds, the quantum discord is con- 
stant and only classical correlation is lost. When f > f', the 
classical coiTelation does not change with time and the quan- 
tum discord decreases. 

Fig. 7 show the transition time is relevant to a. One can 
observe the transition time displays a parabola against time 
and the symmetry axis locates dXa - a' . The transition time 
increase with a in the interval [1, a'] and decrease in \a' , 1]. 
The three-site interaction can shorten the transition time t' in 
most cases compare with the case of or = 0. In some special 
intervals, the interaction can prolong the transition time. This 
feature can be explained by the analysis employed in Fig. 6. In 
the special interval of a, the decay of quantum discord can be 
delayed, so it is reasonable to expect quantum discord could 
keep constant for a longer time than that with other a. Namely, 
the transition time can be prolonged in the special interval. 
This is important because a stable quantum discord do favor to 
quantum information processing. A. Brodutch et al. llsoll have 
proposed that the altering of quantum discord is an indicator 
of failure to a local operations and classical communications 
in implementation of the quantum gates. So it is Valuable to 
keep quantum discord a stable value which could be realized 
via prolonging the transition time. 

In Fig. 8, we adopt fitting method to show the depen- 
dency relationship between transition time and the degree of 



anisotropy y. we find t will decrease when y goes from zero to 
a critical point around y = 0.6 and increase when y goes from 
the critical point to the Positive axis direction. At the critical 
point, the transition time is reach a minimum value. The plot 
illustrate the transition time can be shorten by the degree of 
anisotropy y. 



IV. CONCLUSIONS 

In summary, we have studied the time evolution of quan- 
tum discord and entanglement for two-qubit coupled to a spin 
chain with three-site interaction in the weak-coupling region. 
We analyzed the quantum discord and entanglement for two- 
qubit prepared in a class of X structure state. If the system 
evolute from pure state, we found there exist certain special 
interval in which the decay of the quantum discord and en- 
tanglement can be delayed. At the critical point of the envi- 
ronment, the quantum discord and entanglement suffer sudden 
death at the same time. If no transverse field is in presence. We 
found the quantum discord and entanglement displays oscilla- 
tory decay with time, the absolute value of amplitude is very 
sensitive to the three-site interaction. We also considered that 
two-qubit initially prepared in Bell diagonal state. The de- 
cay of entanglement is delayed in the special interval as in the 
case of pure state, while the decay of quantum discord is en- 
hanced. We have observed sudden transition between classical 
and quantum decoherence for this case, it is found the transi- 
tion time can be lengthened at the special interval of three-site 
interaction and be shorten by degree of anisotropy. Besides, 
we have considered the eff'ect of anisotropy parameter on the 
decay behaviors of quantum correlations. It has been shown 
the effect of three-site interaction on quantum discord can be 
weaken by the anisotropy parameter. 
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